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Abstract: We use perturbative QCD to calculate the parton level cross section for 
the production of two jets that are far apart in rapidity, subject to a limitation on 
the total transverse momentum Qq in the interjet region. We specifically address 
the question of how to combine the approach which sums all leading logarithms in 
Q/Qo (where Q is the jet transverse momentum) with the BFKL approach, in which 
leading logarithms of the scattering energy are summed. This paper constitutes 
progress towards the simultaneous summation of all important logarithms. Using 
an "all orders" matching, we are able to obtain results for the cross section which 
correctly reproduce the two approaches in the appropriate limits. 
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1. Introduction 

Events with high-pr jets separated by large rapidity gaps at hadron colliders are of 
significant theoretical interest. Studying these events provides us with the possibility 
to better understand QCD in the high energy limit and also to understand QCD 
radiation in "gap" events. There are two major approaches to the production of two 
gap-separated jets. 

In the first approach, parton-parton elastic scattering with a QCD colour singlet 
exchange is regarded as providing the leading contribution to the cross- sect ion. The 
leading- F terms (Y is the rapidity interval between the jets) are summed, i.e. 



(i.i) 



This is the so-called BFKL [Q approach to the gaps between jets process. The ob- 
servable calculated in this approach does not consider any radiation into the interjet 
region. Experiments though, impose an upper bound on this radiation by necessity. 
Calculations 0, [| |], || [J based on this approach have been performed and compared 
with some success to the data from HERA and Tevatron |7]] . 

In the second approach the energy flow into the gap is certainly considered. In 
accord with experimental observations soft radiation with transverse energy below 
Qo is allowed in the interjet region. This gives rise to logarithms of Q/Qo where Q is 
the transverse momentum of the jets. The aim here is to sum all terms of the form 

~ a n s Y m L n , m<n (1.2) 

where L = \nQ 2 /Q 2 ) . In particular, at each order of a s the full F-dependence of 
the coefficient of L n is kept. The global leading logarithms of Q/Qo (LLQ ) have 



been summed for various jet definitions || |9], [T^, [TTJ] and non-global effects have been 
considered in JCT], [T2|] . 

For asymptotically large Y at fixed L the BFKL approach is expected to be the 
most appropriate one. On the other hand, using a gap definition with large values of 
Qo at realistic collider energies the second approach ought to be the best. In order 
to get a better understanding of the gaps-between-jets processes at colliders it is 
desirable to combine the two approaches. This is the main issue of this paper. 

To combine the two approaches order-by-order we need to prevent double count- 
ing and make sure the divergences arising from the BFKL approach (at each order 
in a s ) cancel in the jet cross section. To this end we recalculate the gap cross section 
with limited transverse energy-flow into the gap in the high energy limit. In addition 
to the double-leading-logarithmic (DLL) terms we also consider certain terms sub- 
leading in Y, namely those arising from the imaginary parts of the loop integrals. In 
the first instance, we thus sum terms of the type 

~ a^L n Y rn n n - m , m < n (1.3) 

where only integer powers of n 2 appear. We term this approximation the n 2 DLLA . 
We argue that the gap cross section can be obtained by calculating the elastic quark- 
quark amplitude with all virtual gluon rapidities constrained to the gap region and 
all transverse momenta constrained to be above Qo- Based on this we show that 
in the n 2 DLLA , and also in the LLQ A, the BFKL cross section can be included 
explicitly up to 0(oP s ) and derive the expression for the combined gap cross section. 

At this stage, we have not been able to go beyond 0(a b s ) in fixed-order matching. 
However, we do present an all-orders cross section that smoothly interpolates between 
the LLQq and the BFKL results. In this context, we consider three phenomeno logical 
matching schemes. 

The remainder of the paper is structured as follows. In the next section we 
calculate the gap cross section in the n 2 DLLA , consider some approximations to 
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it and compare to the LLQq cross section. Then, we calculate the quark-quark 
scattering amplitude in the leading-!" (BFKL) approximation, deriving the general 
result at any fixed order in a s . Finally we discuss the matching of the two approaches. 
Firstly we attempt a fixed-order matching which we are able to accomplish to 0{a b s ) 
and then we present our all-order matching results. 



2. Summing logarithms in Qq 

We seek to calculate the cross section for two-jet production in the high-energy 
(i.e. high rapidity separation) limit, with limited total scalar transverse momentum 
in the interjet region. We require this transverse momentum to be below Qq, the 
jets to have transverse momentum Q, and be separated by a rapidity interval Y, and 
consider the region 

^qcd « Ql « Q 2 « s = e Y Q 2 . (2.1) 

The first of these relations is to ensure that we can use perturbation theory, although 
we will ultimately be interested in extrapolating Q towards Aq C d, while the other 
two are to ensure that the logarithms we are interested in are large. In this section, 
we calculate the partonic gap cross section for the process, qq' — > qq' 

° = dQ 5 ( 2 - 2 ) 

for fixed s in the 7r 2 DLLA . Since we are not sensitive to collinear emission, we do 
not need to worry about the details of initial state collinear factorization, and only 
calculate partonic cross sections. 

Our approximation implies the eikonal (soft gluon) approximation, valid when 
Q 2 <C s. To generate the leading logs in Q , we make the approximation of strongly- 
ordered transverse momenta for the real and virtual gluons. The scalar sum of 
all gluon transverse momenta can therefore be approximated by the largest single 
transverse momentum, ^ q±i maxj q±i. We therefore require that any real gluon 
emitted in the gap region have transverse momentum below Q . 

Let us denote denote by a^> the gap cross section at 0(a™ +2 ). Clearly at zeroth 
order, the gap cross section coincides with the inclusive cross section, 

,(„) _ ^ (N? - 1) 

As the basis for the calculation to all orders we employ the following theorem. 
2.1 Theorem 

"At any finite order in perturbation theory, in the strongly-ordered approximation, 
the gap cross section is given by the two-to-two cross section, with all virtual gluon 
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rapidities constrained to the gap region and transverse momenta above Qq, and with 
the exchanged gluons attached only to the external lines." 



This is clearly a major simplification, since it means we never have to calculate 
any real emission or triple-gluon-vertex diagrams. 

Central to the proof of this theorem is the following result. For every cut through 
a given cut diagram, and every way of attaching a softer gluon to the external lines 
of that diagram, we can always make four cuts: in which the soft gluon is to the left 
of the cut; to the right of the cut; crossing the cut from left to right; and crossing 
the cut from right to left. The contributions from the real parts of the loop integral 
are equal and opposite to the real emission integral. The imaginary part of the loop 
integral, if present, cancels between the cuts to the left and right of the gluon. The 
total contributions are therefore equal and opposite. 

It is then straightforward to see that we only get a non-zero contribution from 
configurations in which all gluons, real or virtual, have transverse momentum above 
Qq. If the softest gluon is below Q , then the value of the observable is the same, 
whether it is real or virtual. The cancellation is therefore able to take place and we 
get zero contribution. Therefore the softest gluon must be above Qq, therefore all 
gluons must be above Qq. 

We can also use the same argument to prove that the softest gluon must be 
within the gap region. The same argument does not however apply to harder gluons, 
as it is not guaranteed that for them all four cuts listed above are consistent with 
strong ordering. We assume here, as in conventional calculations of gaps-between-jets 
processes, that contributions from outside the gap region cancel, but note that our 
recent work has cast doubt on this assumption [13]. Therefore, all gluons must be 
within the gap region. Finally, since for a real gluon in the gap region our observable 
is zero, all gluons must be virtual. 



2.2 Calculation in the n 2 DLLA 

We therefore have to calculate the all-orders amplitude for quark 2 — > 2 scattering, 
with the phase space for the gluons constrained to the gap region in rapidity and 
with transverse momentum above Qq. The loop integrals themselves are trivial - 
they are just nested sets of integrals. The only complication concerns the colour 
structure. 

Let us denote by AC(Q ) the amplitude defined above, including its colour 
factor. It can be decomposed into singlet and octet components, 

AC(Q ) = AtiQo) C a + As(Q ) C 8 . (2.4) 
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The gap cross section we are interested in is therefore proportional to 



a oc |A(go)| 2 Tr{C t 1 C 1 } + |^ 8 (g )| 2 Tr{CtC 8 } = {A^Q^N* + \A 8 (Q )\ 2 ±(N 2 C -1). 

(2-5) 

We denote the nth order contribution to the colour amplitudes as A\ n \ Since for 
Qo = Q the transverse momentum integrals give zero, we must have 

A l {Q)=Af\ (2.6) 

Since the zeroth order contributions Af^ are given by single gluon exchange, we 
have A^ = 0. We choose to normalize ^4g ^ in such a way that the constant of 
proportionality in ( p.5[) is unity, i.e. 



We introduce a simple vector notation, with basis vectors e x = Ci and e 2 = C 8 , 
such that 

A=( A A, (2.8) 



and we can write 
with 



As, 

<t = A\Qo)SA{Qo), (2.9) 
N 2 



o (2 ' 10) 

We explicitly calculate the amplitude at 0(a 2 ) and 0(a^). Note that, besides 
the (real) leading- parts of the amplitudes we also keep the sub-leading terms which 
differ from the leading parts in the powers of Y and tt, see appendix. As for the one- 
loop amplitude A^(Qo) the soft gluon exchange changes the colour of the external 
quarks, so that the one-loop correction can be considered to be a matrix acting on 
,4 (0) to obtain A (1) . Similarly, A {2) can be obtained from A^ 1 '. We can therefore 
write the result in matrix form as 

A^(Q ) = - 2 -^ f Q ^TA^, (2.11) 



71 



and 



A {2) (Qo) = ( 2 -^) 2 f ^r 2 # 



71 



1 f 2a s f Q dk±. N ' 



2! V tt J Qo k A 



A {0) (2.12) 
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with 



r=( ° )• (2-13) 



This can be summed to all orders by inspection, to give 



^(Qo)=exp|-^^ f^j^)' ( 2 - 14 ) 

where the exponential of a matrix is defined by its power series. This result can also 
be derived by considering the evolution equation 

A(Qo) = A(Q) - ^ f ^TA(k ± ), (2.15) 



which can easily be shown to be solved by the expression above. 
We therefore have 

aU ^ ( 2a s f Q dk ± ,1 f 2a s f Q d&, 1 t ,„. 
a = A^{Q)exp{ -/ -^rHSexp^ -/ -^rUQ. 2.16 

As a final step, we note that a scalar can be written as its own trace, then the 
cyclicity of the trace can be used to bring A to the beginning, and finally, we can 
define 

A{Q)A\Q) = H = ( ° a i ) , (2.17) 

V i(^ 2 -i) 



to give 

2a, f Q dk ± A „ f 2a, f Q dk A 



a = Tr 



Hex P (-^r^rt)sexp(-^r^r 

I 7T Jq q k ± J [ 7T fc ± 



(2.18) 



If we neglect the imaginary parts of the anomalous dimension matrix, it becomes 
diagonal, so that the matrix structure collapses and the result is trivial, 

a crW exp {-N C ^LY) ee <j dl , (2.19) 

where L = hi^p. (|2.19|) is the gap cross section in the DLL A since the neglected 
imaginary parts of T are the terms sub-leading in Y. 

However, we are interested in the solution with the imaginary parts. This can 
be obtained by diagonalizing the anomalous dimension matrix T, since then the 
exponentials become straightforward. We define R to be the matrix that diagonalizes 

r, 

fEER-rR=( A ^ ) A ° ) ), (2.20) 
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where A* 1 ' 2 ) are the eigenvalues of T, which we order by 

He(A (1) ) < He(A (2) ). (2.21) 
We can therefore insert factors of RR -1 or R -1 tRt at appropriate points, to obtain 



a = Tr 
= Tr 



R^HR- 1 * exp [-^LRtrtR-^J RtSRexp {-^LR-^r}] (2.22) 



HD f SD 



(2.23) 



where 



H = R^HR -11 (2.24) 
S = R f SR (2.25) 

»-{-?"'}-(-' { r , „ ( -i a - l )' **> 



For the eigenvalues, we obtain 



\(i) _ N c Y - 2N c™ - N W N c Y2 - 4Yin - 4tt 2 
N 2 - 1 vr 2 N? - 1 in 3 ^ ( 1 



(2.27) 



2iV? Y N* Y 2 \Y 3 



A ( 2) = jgy - gA^r + Agyggg - - 4tt 2 (2 ^ 2g) 
2 iv, v*v 

where the approximation gives the leading real and imaginary parts in the large- F 
limit. 

For the diagonalization matrix, we obtain 

R = -4iNl*m ; ■ (2.29) 

It is worth noting that this becomes diagonal in the large- F limit (with leading 
off-diagonal terms purely imaginary and of order 1/Y). 

In terms of A^, the matrices S and H have fairly simple forms, 

s = [ m _ ^ ) • ( 2 - 3 °) 

~ _ 4q(°)(iV c 2 -l)7r 2 / |A«| 2 -mA« 

((N 2 - 1)tt 2 - 4AW 2 iV2) ((iV 2 - 1)tt 2 - 4A«* 2 iV c 2 ) V^ (1) * ^ 

(2.31) 
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121 



112 



The result for the gap cross section reads 

a= exp{-^L(A« + A«*)}s 11 H 1 
+ exp{-^L(A( 2 ) + AW*)}s 12 H 2 
+ exp{-^L(AW + A( 2 )*)}s 21 H 1 
+ exp { -°^L (a< 2 > + A^*) } S 22 H 22 . (2.32) 

Since S and H are Hermitian, it is easy to see that this result is purely real, as it 
should be. We do not explicitly substitute in for A^, as the results do not particularly 
simplify further. 

To obtain the asymptotic large-!" limit (i.e. Y ^> tt), we take the large- Y" limit 
of the exponents and coefficients separately, 

( a s T N? - 1 7T 2 1 N 2 - 1 TT 2 

a = o-(°) exp <! -L 



7T N? Y J iV c 4 Y 2 
a, i N} - 1 TT 2 



_ a (0) e xv \_^ LNcY }4^ 
F l 2vr J 1 



N A Y 2 



+(7 (o) exp J -^ LNcY 1(1 + 3^—^ ) . (2.33) 



N 2 r - 1 TT 2 \ 



tt J V Nt Y 2 J 

Clearly, the final term contains the expected exponentiation of the lowest order term. 
However, for large Y or large L, the first term dominates. Indeed for large enough 
Y, the result is actually Q - and a s -independent, 

2.3 Extending to the LLQ A 

It turns out that the extension of the previous results to the LLQqA, i.e. including 
those additional terms which lie beyond the high energy approximation, is straight- 
forward. The basic structure, ( |2.18|) , remains unchanged, but one has to modify the 
hard and anomalous dimension (but not the soft) matrices. The hard matrix H is 
sensitive to the kinematics and dynamics of the hard scattering, for example, the 
quarks' spin and flavours (whether they are an equal or unequal flavour pair). The 
anomalous dimension matrix is sensitive to the hard scattering kinematics (i.e. the 
scattering angle) and also to the definition of the gap observable, in particular the 
shape of the region over which the scalar transverse momentum is summed. 

In our case, since we consider a process that has only one colour flow at lowest 

order (qq' — > qq') the hard matrix is still proportional to ( ^ j and the difference 
appears only as a change in (see ( |2.17| )) which is an overall factor that multiplies 
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all of our results. In the following we denote this changed LO cross section by a^ 0> 



full' 



The anomalous dimension matrix, T, is a function both of Ay (the rapidity interval 
defining the gap) and Ar] (the rapidity interval between the two jets). In the cone 
algorithm, one has Ay = Arj — 2R, where R is the cone radius used to define jets. In 
this case it is natural to define Ay = Y, and obtain 



VK 



■VK 



iY, 



v 1_ • . 



(2.35) 



where 



p(Y) = log 



log 



sinh(Ar//2 + Ay/2) 
sinh(Ar]/2 - Ay/2) 
sinh(F + R) 
sinh R 



Y 



y Y -^-io g (i-e- 



2_R\ A_largo Q -2R 



(2.36) 
(2.37) 



R = 1 is common for this kind of study (although R = 0.7 has also been used). Note 
that p(Y) tends to zero for small Y and is monotonic, so always remains small. In 
the LLQqA the coupling has to run. This can be easily incorporated into the cross 
sections obtained so far, simply by replacing 



a 



7T 7T 



a s (kx) 



In 



1 



f3 a s (Q)L : 



(2.38) 



where (3q is given by 



(2.39) 



The additional terms in the matrix V are proportional to the identity matrix. 
Therefore, the diagonalization matrix R and hence also H (with a modified expres- 
sions for cr(°\ see above) and S ( |2.30| , |2.31|) are unchanged. The only effect is to add 



a constant N ° 1 p(Y) onto both eigenvalues in the matrix D (|2.26|). This results in an 



overall factor in the gap cross section (see fl2.32j )). The full LLQq gap cross section 
(without non-global effects), <Tf u ii, thus reads: 



Ofuii = o exp 



N 2 - 1 



In 



1 



Poots{Q)L 



(2.40) 



a is the same as a in the previous section but with a running coupling and a 



: , , _ n jo) 

replaced by crf^ u . The non-global logs |14, |T^] have not been resummed in the case 



of our observable, yet. An estimate of the effect for a different gap definition can be 



found in 11 



It is worth noting that the structure of (|2.40|) is true for all flavours of hard 
process. However, it is only true if the gap definition is azimuthally symmetric and 
one obtains a more complicated structure for the gap definition used by HI and ZEUS 



- 9 - 



and calculated by Appleby and Seymour JTT j , or for an rj-(j) 'patch' as considered by 
Berger, Kucs and Sterman |T0| . 

Fig.[l] compares the DLLA (dotted line), the n 2 DLLA (solid line) and the n 2 DLLA 
with running coupling effects included (thick solid line). All results are taken in ratio 
with the full LLQ A except that the lowest order cross-section is always taken to be 
cr(°\ even in Gf u u (i.e. we are not concerned here with the small difference between 
cr(°) and crfuii)- We use R = 1 and a s (Q) = 0.14. The running of the coupling clearly 
plays an important role. 

FigfJ compares o with a running coupling to <Jj u u at different values of L (again 
with the lowest order cross-section taken to be in both cases). The tt 2 DLL result 
increasingly differs from the full LLQ Q result towards larger values of L. One obtains 
the full cross-section to rather better than 10% for L < 7. 



L=4 



a 



full 




Figure 1: The ratio of <Tf u u/a x for R = 1 



2.4 Numerical results 

Before we continue to the BFKL cross section we numerically compare the cross 
sections a, cr DL and in ( |2.32| ), ( |2.19| ) and ( |2.34| ) using a fixed coupling, a s = 0.2. 



FigfJ shows these quantities (divided by a^) for two values of L. a^, obtained by 
taking a in the large-V limit, is divergent at Y — > and independent of L. This 
is in contrast to a which falls with increasing L at low and intermediate values of 
Y . As a result, the very good agreement between a and at quite low values of 
Y, as illustrated in the lower pane of Fig. 3, is somewhat accidental. The agreement 
becomes worse again as L increases. Of course, the two are always coincident for 
large enough Y. The cross section ctdl is strictly leading-y but, as demonstrated 
in ( |2.33| ), it is not the large-F limit of o. Contrary to a^,, the approximation of a 
by cfdl gets better with decreasing L. So, including the imaginary non-leading-F 
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0.91 - '■ 7 

0.9 I 1 1 1 1 

1 2 3 4 5 

Y 



Figure 2: The ratio of o^nja for R = 1 and a running a s 

parts of the loop integrals has completely changed the physical result in the large Y- 
and/or L-limit. 



3. Matching with BFKL 

In this section we derive the elastic quark-quark scattering amplitude via BFKL 
(colour singlet) exchange at any fixed order in a s . Then, we determine the terms 
it has in common with the ir 2 DLL amplitude calculated earlier. We discuss how 
to combine the two to provide a cross section that includes the n 2 DLL gap cross 
section and the first orders of the 2^2 BFKL cross section without double counting. 
Finally, we discuss a more phenomenological way to combine them to all orders to 
provide a cross section that agrees with each in its domain of validity and smoothly 
interpolates the two. 

We work in Feynman gauge and in d = 4 — 2e dimensions. We continue to use 
the normalisation of ^.g '' as given in ( |2.7| ) with a modified expression for <To though. 
In 4 — 2e dimensions <r Q2.3p is replaced by: 

a^-^(^X-^- (3 1) 

^ ~^ \ Q 2 ) r(i-e)" {6A) 

3.1 The BFKL amplitude 

The elastic quark-quark scattering amplitude in the leading-F approximation can be 
written as 

A BFKL (Y q 2 ) = ^sal (2^ / ^n^t /( k i» ^ * Y)> (3-2) 



AT? J k?(ki-q) 
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L=2 




0.01 r 



0.001 



0.1 r 



0.01 r 



0.001 



L=5 




Figure 3: The gap cross sections normalized to at different values of L (a s = 0.2). 



where / is given by the BFKL equation |l 




dY 



/(k x , k 2 , q; Y) = 5 2 - 2e (k! - k 2 )5(Y) + [u(kl) + ^((k! - q) 2 )] /(k 1; k 2 , q; Y) 



(27T)!- 



2c 



Ci 2 " 2 ^ 



q- 



I /k 2 + (kx - q) 2 



£l(£-q) 2 (ki -£) 2 V^ 2 (*-q 



/(^k 2 , q; y) 

(3.3) 



with 



^(k 2 ) 



1 Qf. 



rf 2 ' 2e l k 2 

2(2tt) 1 - 2£ 7 £ 2 (£-ki) 2 



a. 



7T 



(3.4) 
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In these equations we have q 2 = Q 2 . We present two solutions for the elastic ampli- 
tude in the leading- F approximation. First, we solve (|3.3|) by iteration and derive 
the general expression for the amplitude at 0(a™). The second solution has been 
derived in 0]; the result is written as a sum over conformal spins. 

In accordance with our previous notation, we denote by the contribution 
to / proportional to a™ -1 (which therefore contributes to the n + l-th order of the 
amplitude since the quark impact factors each contain a factor a s ). The zeroth order 
expression in ( |3.3j ) is therefore The momentum dependent part of the elastic 
amplitude (CO) reads 



d 2-2e kid 2-2e k: 

k?(k x - q) 



a s Y 



n(2vr) 1 - 2e 
1 



V (n) (ki,k 2 ,q ; y) 

d 2-2c ki d 2-2e k2 d 2-2e £ 



+ 



2(£-ki) 2 V^ 2 (ki-q) 2 k 2 (£-q) 2 

,2 



/(^(kj.k^qjy) 



+ 



q- 



L£ 2 (£-q) 2 k 2 (k 1 -q) 
1 / 1 



+ 



1 



(^-kx) 2 V^ 2 (ki -q) 2 k 2 (£-qf 

! k ' ^(ki,q;y) 



/ M (^k 2) q;y) 



k?(k!-q) 



(3.5) 
(3.6) 



In the second term of (|3.5|) we substitute £ <-> kj.. The factor that multiplies /( n_1 ) is 
invariant under this transformation. The only effect therefore is that the arguments 
of the functions /( n ~ 1 ) in the first and in the second term become equal and we can 
combine the two terms. Note that this is only possible since the quark impact factors 
are constants. Performing the k 2 integration we end up with the following simple 
resursion formula for the function as defined in ( |3.6| ) 



ff w (k liq; y) 



]2-2e. 



n(27r) 1 - 2e 
2 1 



q 



_£ 2 (£-q) 



kf 

2 (£ - k x ) 2 V £ 2 



2n 



x (Ait) 



r(i + e)r(- 

T(-2e) 



,2\-€ 



2\-e 



(ki - q) ; 



;((ki-q 



g^- 1 \k h q;Y) 



,2\-£ 



3 M (ki,q;r) 



(3.7) 
(3.8) 
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To obtain the amplitude at we have to take the bracket in (3/7) to the power 

of n— 1 which gives rise to a double sum with two binominial coefficients and we have 
to perform the \z\ integration. Our result for the elastic quark scattering amplitude 
is 

.(») _,7 iV»>l<°)„ ( as Y~ l vn-i N 2 c -l ( q 2 



Hi + ^HT -^f/n-nm a(i + aq 

r(-2e) ) ^{ k ){j) { - 1/2) A(7)A(*-i) ( } 
with 

iYl + xe) /n\ n! 



1 + x)e) ' V^/ fe!(n- fc)!' 



A(X) = 

To simplify further we perform the summation in ( |3.9|) for the first the few orders 
in a s . The result can be expressed in terms of functions j n (e): 

*w-SS (3 ' 10) 

which obey also the following relation: 

Note that 7 n (e) diverges as 1/e. Suppressing the argument of the functions 7 n (e) we 
obtain: 

A^-^^f^);. ~i (3.12) 



.(2) _ , ,(0) gg AT c 2 -l /_q^ x " 

(2tt) 8iV c Uvr/x 2 ; ,uwu ,iy e 2 



^L = ^7^T-hi- y iTTra) 7o(7o-7i) ~ -, (3.13) 



4( 3 ) _ , ,(0) «| ^c-i y2 f_q? x !i 

7o ryo(7o - 7i) - 7i (7o - 72) + ^72(70 - 7i)J ~ ^ ( 3 -14) 

Here, we have also noted the leading behaviour of the amplitudes as e — > 0. 

The all orders leading-F result for the elastic quark-quark amplitude, Abfkl, 
is written as a power series in a s in our approach. In contrast to this it can also be 
written as a sum over conformal spins @] 



Abfkl — — iA s 



c 



2tt A" 2 



00 

X 

m=— 00 



£ / ^ ^ p + (m - 1/2)^ + (m + 1/2) 2 ] eXP [ " 2m(;/) F] ^ (3 ' 15) 
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with 



uJv) = ^-^ [2^(1) - V(l/2 + M/ 2 + iz/ ) - ^(1/2 + ln.l/2 - iv)\ (3.16) 

7T 

and ip{z) = 1/z dT(z)/dz. It is worth emphasizing that this all-orders expression 
is finite, even though its expansion to any finite order is divergent, as seen above. 
Physically, this is because the divergent terms exponentiate to zero. This property 
will be crucial in constructing an all-orders matched cross section below. 
In the limit Y — * the BFKL 2^2 cross section reduces to [|J 

N 2 — 1 7l 2 

(?bfkl\y-+o = N c \-Abfkl\ Y ->o = cr( ° ) V 4 Y 2 ' 

This behaviour will also be important in the all-orders combination of BFKL with 
the gap cross section. 

3.2 Singlet exchange in the ir 2 DLLA 

As the first step towards a matching of our gap cross section to BFKL we extract 
from the elastic amplitude in the ir 2 DLLA its singlet exchange part taken in the 
leading- F approximation. To be more precise, we extract from the resummed elastic 
amplitude (|2.14j ) the singlet exchange contribution A^iQo) at fixed order in a s . 



Keeping in this expression only the terms leading in Y we obtain an expression we 
term Aig(Qo). The subscript stands for 'strongly ordered' since at Qo = this 
amplitude differs from the leading- F (BFKL) amplitude in that it is calculated in 
the strongly ordered approximation. For the first few values of n we have: 

(A {0 \Qo)Ks = o, (A {1) h,s = i* ^4 0) i--L) - a (— — L) 

4iV ( f 7T 7T 

U {2) (Qo)hs = h-^L) 2 (T 2 A^\ = \a (— — L) 2 (^Y), 

21 7T 2 71 2 

(A {3 \Qo)Ks = k-^Lf(T^% = U-^Lf {^Yf 
3! it 3! 7r 2 



The general result for n > 1 can be written as 

and one obtains for the resummed (LY) singlet exchange amplitude 



A 1 ,s[Qo) - -« 2iy3 y 8 



c 



oxp ( YL 



(3.19) 
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In the limit of large Y or L the exponential vanishes and we obtain a finite result for 
the all-orders amplitude Ai } s(0)- 

A hS (o) = -4i,s(Qo)|_ = -i y A * ] - (3 ' 20) 



Using (|2.7|) the corresponding cross section as reads 

o-(o) 2 2 

ClTXpCl 1^(0)1 = - m y 2 



a s = C\ \A^)\ 2 = t^C? |A, 5 (0)| 2 = a® %^ (3.21) 



This is exactly the large Y limit of our gap cross section ( |2.34| ) and the small Y limit 
of the BFKL cross section ( ]3 . 1 7] ) 



&bfkl\y^o — °~s — °"ly^oo (3.22) 

We will exploit this remarkable fact to construct an all orders combined cross section 
further below. 

3.3 Fixed order matching to BFKL 

In order to combine the 2^2 BFKL cross section with the gap cross section a we 
have to ensure that a) at each order in a s the e-poles in the BFKL cross section are 
cancelled and b) we have avoided any double counting. The only contribution to a 
which also appears in the BFKL cross section is -4[ n ](0) (J3TTJ). Let us keep Q ^ 
for the moment and consider the generalization of the amplitude (|3.18| ) to d = 4 — 2e 
dimensions: 



2iV c 3 Y n\ \ 7T r(l - e) J Qo k 1+2 ' 

We have calculated this amplitude explicitly in d dimensions only up to 0(a^) accu- 
racy (see appendix) , but using the well-known fact that the eikonal Feynman rules are 
unchanged in d dimensions, so only the phase space integrals become d- dimensional, 
it is straightforward to see that ( |3.23| ) will hold to any order. Now we can set Qo = 
and obtain the elastic scattering amplitude in the leading- F approximation and in the 
approximation of strongly ordered transverse momenta (indicated by the subscript 
'5'): 



Now, we calculate the difference of the BFKL amplitude ( |3.9|) and its strongly ordered 
approximation ( |3.24|) for the first orders in a s : 



5 {n) = AIkl ~ 41(0) (3.25) 
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^) = -^4V^lc(3)6 2 + 0( e 3 ) (3.26) 

6(2) = iA ^l ¥ ^ C(3)e + ° (e2) (3 ' 27) 

5® = -iAfa^ ^ (iV c 2 - 1)C(3) + 0(e) (3.28) 

5 (4) = C(l/e) (3.29) 

Here, we have set e = in A^ . ((x) is the Riemann (^-function. The leading pole 
in the BFKL amplitude is captured by the strong ordering approximation. And, 
not surprisingly, the strongly ordering approximation cannot provide all e-poles; at 
0(al) the difference of the two expressions starts to become divergent. Working in 
the ir 2 DLLA (or in the LLQ A) will allow us therefore only to include the BFKL 
cross section to a fixed order of a s . Note that the smallest power of e in Ag FKL 
is e _1 , in A§ FKL it is e -2 etc. (see ( p.l2|) -( j3~T4| )). Remarkably, the strong ordering 
approximation therefore agrees with the full amplitude not only at the smallest power 
of e but rather at the three smallest powers. 

Now we come to our original goal to include the BFKL cross section in our gap 
cross section. We start from the Theorem stated above and work at fixed order in 
a s . Denoting the production amplitude (in the n 2 DLLA , including color factor) for 
more than 2 particles by B(Qq) the theorem reads 

= |A(0)| 2 C 2 + |A(0)| 2 C 2 + \B(Q )\ 2 = |A(Qo)| 2 C 2 + | A S (Q )\ 2 C 2 8 (3.30) 

where the squares are to be read symbolically representing the sums over A^^A^ 
(and B(Qo) , respectively). In this paper, we wish to also include the colour singlet 
BFKL 2^2 amplitude A^ FKL . 1 However, A\(0) also includes terms subleading 
in Y , whereas the BFKL amplitude is leading-!^. We therefore have to keep these 
subleading terms; they are given by (Ai(0) — Ai,s(0)) since Al^O) is calculated in 
the leading- F approximation. Omitting indices for the moment we therefore define 
the following fixed order gap cross section 

= \Abfkl + A(0) - A,s(0)| 2 C 2 + | A(0)| 2 C 2 + \B(Q )\ 2 
= \A 1 (0)+5\ 2 C 2 1 + \A 8 (0)\ 2 Cl + \B(Q )\ 2 

= |A(0)| 2 C 2 + I A(0)| 2 C 2 + \B(Q )\ 2 + (2Mb [Ai(0)<T] + |5| 2 )C 2 
= a« + (21Re [Ai(0)5*] + |5| 2 )C 2 

where, in the last line we have invoked the theorem ( |3.3U ). Exploiting the fact that 
5 ( 3.25 ) is purely imaginary and showing the indices again we end up with 

af] p = a {k) + A (fc) , A (fc) = N 2 ^ [2ImAj m) (0) ■ (-^ (n) ) + <f( m W n >] (3.31) 

m+n=k 

1 We do not aim here to include those BFKL logarithms generated as a result of real emissions. 
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It now is easy to see at which order of a s the function A( fc ) and hence a^p starts 
to be divergent. 5^ is finite for n < 3, see (|3.26| )- (|3.28|) . Em^l^^O) only has a 
1/e-pole; this gets multiplied by a constant, namely 5^ Q3.28| ), for the first time in 
A 1 - 4 -*. Hence, A^ k ' is expected to include e-poles for k > 4. To calculate A^ 2 \ A*- 3 -* we 
need A^(0), ^ 2) (0) which are given in the appendix (A^ is not needed since IS 
zero). Comparing with ( |3.24j ) it turns out that these amplitudes have no subleading 
(in Y) parts, 



Em4 fc) (0) = InU^(O) , k = 1,2. 



(*), 



We can therefore use the leading-!" expression 

instead of ^P(O) in ( CT ) 
and inserting the definition of 5^ n ' we arrive at: 



i Mi 



m+n=k 



k < 3 



Nl J2 fwg^W^-Wg(0)Im45(0)l (3.32) 



n+m=k 



Using (pT2|), (|3T3|) and (pi) in (|3~3l we obtain: 

A< 2 > = 0(e) 

A< 4 > = 0(l/e) 



(3.33) 
(3.34) 
(3.35) 



Since A^ fc ^ diverges for k > A our final gap cross section includes the n 2 DLL result 
and the BFKL cross section up to 0(a^), 



k=0 

= o + A^ 



(3.36) 



where the resummed cross section a is given by ( 2.32Q . 

/±( 2 ) corresponds to the squared 1-loop amplitude; the fact that A^ 2 ^ is zero means 
that the gap cross section a fully includes the leading order BFKL cross section (a 
also includes the terms that cancel the corresponding divergence). 

Note that we could include in a terms that are suppressed in the high energy limit 
(including non-global logarithms). This would not affect the value of A^ needed for 
a combined LLQo+ BFKL cross section. Therefore, we can use af u ii, ( 2.40|) , instead 
of a in ( p.36|) and obtain the cross section that combines the conventional (partonic) 
gap cross section Q with the first orders of the BFKL result (not taking into 
account the non-global logarithms). 
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Fig.|] compares the combined gap cross section a gap (|3.36| ) with a. We use a fixed 
coupling a s = 0.2. Note that part of the BFKL cross section is included in a. Fig|| 
therefore shows the difference between the full (up to 0(oP s )) and the partial inclusion 
of BFKL. The importance of BFKL grows towards larger values of Y and/or L since 
a decreases whereas the BFKL cross section is independent of L and grows with Y. 



"gap ® 

a 




Figure 4: The fixed-order combined gap cross section compared to a (a s 
malised to a^). 



0.2) (nor- 



The inclusion of the first orders of the BFKL cross section in the gap cross section 
is found to change the result by a factor ~ 2 for Y = 4 and L = 5. At the Tevatron 
in Run II jets with pt ~ 30 GeV and rapidity separation of Ar] = 6 (i.e. Y = 4) can 
be expected. For Q = 2 GeV such jets corresponds to L = 5.4. The impact of the 
singlet exchange cross section might therefore be observable at the Tevatron. At the 
LHC values of Y = 6 and L > 5 are most probably accessible, enhancing the effect 
of BFKL even more. 

Fig.|5| shows our results together with the resummed BFKL cross section |4[] . The 
BFKL cross section should be appropriate for large enough Y whilst (Jf u u is the most 
reliable for small values of Y . A combined cross section should therefore smoothly 
interpolate between the large-F BFKL cross section and a at small Y. 

3.4 All orders matching to BFKL 

In the previous section we showed that an order-by-order matching of the LLQ 
and BFKL results can only work for the first few orders, because beyond the fourth 
order, the sub-leading poles in the amplitude do not match and the combination is 
divergent. We now show that, using the fact that the all-orders results for Abfkl and 
^1,5(0) are finite, and that they are identical at small Y, it is possible to construct an 
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L=5 



0.1 r 



0.01 




Figure 5: The fixed-order combined gap cross section compared to the full BFKL result 
(a s = 0.2) 



all-orders cross section that does smoothly interpolate the LLQ and BFKL results, 
agreeing with each in its region of validity. In fact we can construct several different 
matched cross sections that all fulfil these requirements while differing slightly in their 
predictions in the intermediate region. We show numerical results from all three as a 
measure of the uncertainty inherent in the matching procedure. We emphasize that 
these are all phenomenological approaches that cannot be proved order by order, 
since they are based on all-orders summations of divergent expressions. 

3.4.1 Simple Matching 

In the previous section, we showed that the combined gap cross section could be 
written as 

(3.37) 



<y 9ap = <y + N 2 c [2HmA(0) ■ {-i8) + \5\ 2 ] 



where we have dropped the superscripts, implying that we are summing each ex- 
pression to all orders. Although all the expressions appearing in the square brackets 
in (|3.37|) are divergent at every order of perturbation theory, they all have finite 
all-orders expressions. In particular, ^4i(0), the LLQ singlet amplitude for fixed Y 
and L —>■ oo, is zero. We therefore have 



(T g a P = cr + N 2 \5\ 2 . 



(3.38) 



This expression achieves our goal of having a smooth matching of the two all-orders 
cross sections, in that for small and large Y it agrees with the LLQq and BFKL cross 
sections respectively. For small Y (and any L, since they are both L- independent), 
we already showed that Abfkl and Al,,s(0) become identical. Therefore 5 becomes 
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zero and a gap becomes equal to a. For large Y and any L, a and A,s(0) both become 
zero and hence a gap becomes equal to the BFKL cross section. 

3.4.2 Cross Section Matching 

In the previous section, we also showed that for the orders for which it is finite, ( 3.37|) 
is identical to 

(T g a P = ° + N 2 c [ImA B FKL^mA B FKL - XmAi,s(0) Wli,s(0)] (3.39) 
= a + a bfkl ~ cr s - (3.40) 

We take this as the definition of another matching scheme, the cross section matching 
scheme. It has a very simple interpretation: the matched cross section is equal to the 
sum of the LLQ and BFKL cross sections, with the double-counted terms subtracted 
off. Note that according to the discussion in the previous section, this can be written 
in an identical form to (ft .371) , but with Ai(0) replaced by Ai } s(0)- Since At becomes 
equal to Ai } s for sufficiently large Y, and the 5 term that these terms multiply 
becomes zero for small Y, these expressions are equally good within their regions of 
applicability. 

This expression therefore also achieves our goal of having a smooth matching 
of the two all-orders cross sections, in that for small and large Y it agrees with the 
LLQq and BFKL cross sections respectively. 

3.4.3 Amplitude Matching 

Finally, inspired by the form of the cross section matching scheme, we consider a 
similar matching, but at the amplitude level, 

° 9ap = \(N? - l)\A 8 (Qo)\ 2 + N 2 c \Ai, 9 a P (Qo)\ 2 , (3.41) 
Ai >gap {Q ) = A 1 (Q ) + 5 (3.42) 
= Ax(Q ) + A BFKL - A,s(0). (3.43) 

That is, we write the singlet amplitude as the sum of the LLQq and BFKL singlet 
amplitudes, with the double-counted terms subtracted off. Expanding this expression 
and rewriting it in terms of 5, we obtain 

o 9ap = o + N 2 C [21mA (Q ) • H<5) + \S\ 2 } ■ (3.44) 

That is, an identical expression to (|3.37f) but with Ai(0) replaced by A\{Qq). Since 
Al(Qo) becomes Qo-independent for sufficiently large Y, and the 5 term that these 
terms multiply becomes zero for small Y, these expressions are equally good within 
their regions of applicability. 

This expression therefore also achieves our goal of having a smooth matching 
of the two all-orders cross sections, in that for small and large Y it agrees with the 
LLQq and BFKL cross sections respectively. 
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3.5 Numerical Results 

We show numerical results of all three schemes in Figs. [], |7|, [| for L = 2, 5 and 6. 
We see that indeed they all achieve the goal of matching the two cross sections in 
the small and large Y limits and the differences are not large in between. The very 
close agreement that one sees at L = 5 is somewhat coincidental. 

We stress again that all three approaches are formally equally valid and equally 
phenomenological, in the sense that they cannot be justified by a finite order-by-order 
expansion. Only with a better understanding of the real-emission contributions in 
the high energy limit, which also have finite all-orders expressions but cancel these 
divergences order by order, can a more formally-justified matching be set up. We 
reserve this for future work. 

4. Conclusion 

Working in the high energy limit we have calculated the (partonic) cross section for 
the production of two jets distant in rapidity and with limited transverse energy 
flow into the region between the jets. Besides the DLL terms, we have summed 
terms sub-leading in Y stemming from the imaginary parts of the loop integrals. 
This allowed us to identify those contributions that are also included in the fixed 
order 2 — > 2 BFKL cross section. We have shown that the leading order BFKL 
cross section, corresponding to the squared 1-loop amplitude, is entirely contained 
in the ir 2 DLLA. We derived an expression that consistently includes the terms of 
the LLQ series and the BFKL series to 0(aP s ) accuracy without double counting. 
In the LLQqA, the inclusion of higher orders of the BFKL cross section in this way 
is not possible since it implies a divergent cross section. 

We have also studied several phenomenological "all order" matching schemes 
that effectively interpolate between the LLQ and BFKL results. Although they all 
yield similar results, the differences between them cannot be resolved without further 
work, specifically understanding the role of real-emission contributions in the high 
energy limit. 

In this paper we have made a first step towards the unification of the two main 
approaches to the "jet-gap-jet" process. The question remains as to how they can 
be fully combined. 

A. Elastic Scattering Amplitudes in d Dimensions 

We calculate the elastic quark-quark amplitude in the high energy limit in the ap- 
proximation of strongly-ordered transverse momenta of the t-channel gluons. We 
work in the eikonal effective theory. The transverse momentum of the softest gluon 
is taken to be larger than Q . Since we need the singlet exchange amplitude for 
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Qo = we consider d = 4 — 2e dimensions. We work in Feynman gauge. Denoting 
the color indices of the initial and final state quark pair by ik and jl, respectively, 
the colour operators for singlet and octet exchange are 

Ci = Sij5 k i , C 8 = -(SuSjk — —SijSki). 

The tree level amplitude reads: 

2s 

Q 

A.l The one loop amplitude 



^ (0) C 8 = - 9 y^ 2 Cs (A.l) 



There are two Feynman diagrams, the box and the crossed box. In each of them 
there are two configurations of strong ordering possible giving the same result, so we 
collect them together: 



AuC n = 2A<® J ^/(fc)e(Q 2 - k 2 )9(k 2 - Ql)C n , (A.2) 

/(*) = i 9 y e {pik _ ie) ^ k + ie)(k 2 + ie y ( A - 3 ) 

AC X = 2A<® J -0^f(k)e(Q 2 - k 2 )0(k 2 - g 2 )C x , (A.4) 

f(k) = -ig\^ {pik _ te){ ^ k P _ 2 te){k2 + te y (A.5) 

We have again extracted the colour factors explicitly, Co and Cx. k is the usual 
transverse part of k (k 2 > 0). The upper limit on k 2 as imposed by the © function 
is a consequence of the strong ordering condition. Simple colour algebra gives 

c °= -jr c » + 4v^ Cl - (A - 6) 

M^)c 8 + ^ic, 

Putting together the terms and performing some simplification, we arrive at 

A (1) C = AuCn+iAxCx (A.8) 



with 



T = 2 / w^ f{k) 0(g2 - k2)0(k2 " g ° } ' (A - 10) 

/Ad k 
j^y[f(k) + /(*)] 0(Q 2 - k 2 )0(k 2 - Ql) (A.ll) 
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We introduce Sudakov variables and perform the integration with respect to two of 
them. The angular integration in the transverse momentum integral is trivial and as 
the result for T and Q we obtain [d — 4 — 2e): 

2a s (47r/i 2 ^ 



7i 
7T 



r(i-e) 

(47T/X 2 ) 6 

ra-e) 



Y 



dk 



U7T 



k l+2e 

Q dk 



Qo 



k l+2e 



(A.12) 

(A.13) 



Note that the real parts from / and / in Q ( |A.11| ) and hence the terms proportional 
to Y have cancelled. 



A. 2 The two loop amplitude 

According to the colour structure there are four groups of Feynman diagrams (Fig.|D|). 
For each group there are four configurations of strong ordering giving the same result. 
In terms of the functions / and /, (|A.5| , |A.3|) the amplitudes read: 



^1 2 )C T = 4^°) 



d d k x f d d k 



(2vr) d J (2tt) 



a x (h, k 2 ) 6(Q 2 - k 2 )6(k 2 - k 2 )6(k 2 - Ql) C a 



x 



a, b, c, d 



with 



a a (h, k 2 ) 
a b (h, k 2 ) 
a c (ki, k 2 ) 
ad{k\,k 2 ) 



f(ki)f(k 2 ) 
f(ki)f(k 2 ) 
Kh)f{k 2 ) 
f(h)f(k 2 ) 



(A. 14) 
(A.15) 
(A.16) 
(A.17) 



We can express each amplitude in terms of T and Q since the B-functions only 
affect the integration over the absolute value of the transverse momentum which 
we left undone in T and Q, ( |A.12|]A~T3D . The decomposition of the sum of the four 
amplitudes in terms of C§ and Ci is straightforward. We finally find for the two-loop 
amplitude (C being the corresponding colour factor): 



He(^C) = A {Q) - 



\m{A^C)=A^ l -\Q\T 



X c 1 + (^ 



Nl + 3 



c, 



(A.18) 
(A.19) 



Note that the sub-leading terms in the real part of stem from keeping the sub- 
leading imaginary part of the one loop result. 
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Figure 6: The gap cross section in the three matching schemes for L = 2 (a s = 0.2) 
compared to obfkl (dots) and a (double-dots) 
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Figure 9: The 2-loop amplitude 
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